Dissipation effect in the double-well Bose-Einstein Condensate 
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Dynamics of the double-well Bose-Einstein condensate subject to energy dissipation is studied 
by solving a reduced one-dimensional time-dependent Gross-Pitaevskii equation numerically. We 
first reproduce the phase space diagram of the system without dissipation systematically, and then 
calculate evolutionary trajectories of dissipated systems. It is clearly shown that the dissipation can 
drive the system to evolve gradually from the 7r-mode quantum macroscopic self-trapping state, a 
state with relatively higher energy, to the lowest energy stationary state in which particles distribute 
equally in the two wells. The average phase and phase distribution in each well are discussed as 
well. We show that the phase distribution varies slowly in each well but may exhibit abrupt changes 
near the barrier. This sudden change occurs at the minimum position in particle density profile. We 
also note that the average phase in each well varies much faster with time than the phase difference 
between two wells. 

PACS numbers: 03.75.Lm, 05.30.Jp, 05.30.Rt 
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I. INTRODUCTION 

The experimental realization of atomic Bose-Einstein 
condensation has stimulated great interest in quantum 
nature of the condensate of weakly interacting bosons. 
Since the condensate is in the macroscopic quantum 
state, a variety of quantum phenomena show up when 
two condensates are linked. If the two condensates are 
separated by a barrier through which atoms can tunnel, 
a bosonic Josephson junction (BJJ) is built and allows 
for the detection of Josephson-like current-phase effects, 
analogous to two weakly coupled superconductors. This 
phenomenon was initially predicted by Javanainen [l[ 
and was further demonstrated by a number of researchers 
using various techniques [34l^|. In addition, the dc and 
ac Josephson effects 043 as well as the Shapiro effect [5[ 
in BJJs have been discussed. 

Furthermore, many theoretical works reveal that the 
BJJ can exhibit new phenomena that are not accessible 
with superconductor Josephson junctions (SJJs). The 
mean-field theory based on two-mode approximation pre- 
dicted the 7r-phase oscillations and macroscopic quantum 
self-trapping (MQST) [111]. The MQST is a fascinating 
nonlinear phenomenon which arises from interactions be- 
tween atoms. Quantum corrections to the mean-field re- 
sults were also intensively discussed [6rJl2|. providing a 
number of fascinating observations, such as collapse and 
revival of quantum oscillations [6j, occupancy of phase- 
squeezing state 0, disappearance of coherence H,lij], and 
destruction of the self-trapped state [lOj . 

In above studies the double- well condensate is treated 
as an isolated system, but actually it is coupled to cer- 
tain thermal cloud and subject to environment distor- 
tions 13 16]. Ruostekoski and Walls studied the energy 
dissipation in BJJs via stochastic simulations in a two- 



mode approximation and showed that the MQST decays 
away [13j |. Wimberger et al. discussed dynamics of the 
double- well condensate subject to phase noise and parti- 
cle loss and indicated that coherence might be increased 
to a maximum for proper dissipation [15]. They argued 
that this phenomenon can be understood as a stochastic 
resonance of the many-particle system. 

On experimental side, the MIT team carried out the 
first double- well experiment and gained interference pat- 
terns proving the coherence of BEC [l?} ■ In 2005, Albiez 
et al. reported the direct observation of the quantum 
tunneling and macroscopic self-trapping for a single BJJ 
[l8j . They measured the density distribution in situ and 
deduced the relative phase versus time from interference 
fringes. A summary of the observed static, thermal and 
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dynamical properties was given in Ref. 
the dc Josephson effect was also achieved 
bifurcation was studied at the transition from the Rabi 
to Josephson oscillations [2l|. 

In this paper, we concentrate on dynamical behav- 
iors of the double-well condensate subject to energy dis- 
sipation based on the time-dependent Gross-Pitaevskii 
(TDGP) equation. The TDGP equation provides a very 
good description of dynamical evolution of trapped dilute 
Bose gase s at temperature far below the critical temper- 
ature |22|, [23. The two-mode approximation for the BJJ 
is an extremely simplified version of the TDGP equa- 
tion. Although it can give a qualitatively description of 
the BJJ system @, H M, d , the two-mode model ig- 
nores the spatial distribution of the density and phase. 
Whereas, these information could be obtained from the 
Gross-Pitaevskii (GP) equation [22]. 

This paper is organized as follows. In Sect. [Til the 
TDGP equation is introduced to describe the double- well 
condensate, where a phenomenological term is added to 
account for the energy dissipation. In Sect. Mil we repro- 
duce numerically the phase space diagram of the double- 
well condensate without dissipation, and make a compar- 
ison with the two-mode model results. Sect. II VI discusses 
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the dissipation effect. A summary is given in the last 
section. 



II. THE MODEL 

A. Time-dependent Gross-Pitaevskii (TDGP) 
equation 

The general formalism of three-dimensional (3D) time- 
dependent Gross-Pitaevskii equation for dilute Bose 
gases trapped in the double-well potential at zero tem- 
perature reads 



ih 



dt 



2m 



V(r) + gN\iP(r;t)\ 



^(r; *),(!) 



where ip(r; t) is the macroscopic wave function at position 
r and time t, m is the mass of an atom, g — Air?i 2 a/m 
is the nonlinear constant with a the s-wave scattering 
length, and N is the total atom number. The double- 
well potential V(r) is given by V(y) = ^{u? x r 2 x + uj 2 r 2 + 
w s r a) + VbCOS 2 (Trr z /qo) [11], where the first term de- 
scribes a harmonic trap with trap frequencies u) XtVtZ , and 
the second denotes the barrier potential. Hereinafter we 



set uj x 



Assuming that the dynamical evolution mostly takes 
place in the direction containing the barrier, the 3D GP 
equation can be reduced to a simpler one-dimensional 
(ID) form 24|. Such reduction is reasonable in the case 
of a strong trapping in the other two spatial directions, 
say, in the x-y plane. Then the ID GP equation in z axis 
only with dimensionless quantities is given by |25| 



1 dr 



d 2 



\^(z;t)\ 



where z — r z /l z , r — tuj z /2 and j3 = lj^_/lo z . The dimen- 
sionless parameters are scaled by the characteristic length 
l z = y/h/JrnuJz) and energy huj z correspondingly. The 
ID potential is reduced to V\d{z) = z 2 + 2wt,cos 2 '{n z/q' ), 
with v h = V h /(hui z ) and q' Q = q /l z . gm = gNm/(l z h ) 
is the dimensionless interaction parameter. This reduc- 
tion greatly simplifies the computation. The wave func- 
tion ip(z; t) satisfies the normalization condition, 



dz\4>(z; t) 



1. 



(3) 



The reduced TDGP equation can be solved by a com- 
monly accepted numerical algorithm called Split-Step 
Crank-Nicolson scheme with space and time both dis- 
cretized (25|. In numerical calculation, t is measured 
in co z /2. If lu z — 2tt x 10Hz, r = 1 corresponds to 
t = 31.8ms. 



B. The TDGP equation with dissipation 

Choi et al. has described a phenomenological damping 
theory based on the TDGP equation 26] . A phenomeno- 



logical parameter 7 is introduced in the GP equation to 
account for the dissipation rate. The reformed TDGP 
equation is written as 



(1-7) 



dr 



d 2 



9idP 



|^(z;t)| 5 



^(z;t). (4) 



The phenomenological parameter 7 characterizes the 
strength of dissipation and its value may be determined 
experimentally. Choi et al. got that 7 = 0.03 by fit- 
ting theoretical results with the MIT experiments dates. 
This approach is comparable to the generalized finite- 
temperature GP equation derived by Zaremba, Nikuni 
and Griffin [H, 



1 dt 



h 2 5 

2m 



+ V(r) + gn + 2gh - tT]*(r; i),(5) 



where n(r;i) = \^(r;t)\ 2 is the condensate density and 
n(r; t) is the noncondensate density. Dissipation arises 
from collisions between the condensate and nonconden- 
sate atoms. 

Ueda et al. have successfully employed this phe- 
nomenological approach to study the detailed dynamics 
of a sudden rotated BEC in a trap with 7 = 0.03 and fig- 
ured out the vortex lattice formation process 28]. This 
approach is also used to study the Bose-Einstein droplet 
in free space (29| , vortex formation in Bose-Einstein con- 
densates in a rotating double-well potential [30( and the 
decay of dark solitons in harmonically trapped, partially 
condensed Bose gases in the presence of phase fluctua- 
tions or dynamical noise 31] . These works prove that the 
phenomenological damping description does indeed pro- 
vide an efficient numerical tool for the research of BECs 
in an open system. 



C. The phase difference and phase distribution 

As discussed in the introduction, we consider two 
Bose-Einstein condensates with repulsive interaction in a 
Quasi-ID symmetric double-well potential, which could 
be experimentally realized by splitting a single well into 
two parts: one cigar-shaped atomic cloud is cut into two 
separated aligned cigars. It is necessary to construct an 
approximate initial state of the double-well condensate 
at first. We set 4>l for the wave packet mostly localized 
in the left well and 4>r the right one. The original trial 
wave function is forged as a superposition oi 4>l and (fin, 



ip(z; t) = ^l(t)0l(z) + iPr(t)(/)r(z), 



(6) 



where tp L (R){r) = y/n£^rje iet -^ T \ n i(R )(r) corre- 
sponds to fraction of atoms on the left (right) side of 
the trap, and therefor n^(r) + tir(t) = 1. They are 
constructed as linear combinations of the two lowest sta- 
tionary modes of the ID GP equation, the symmetric 
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4>+ and antisymmetric </>_. 4>l(z) = ^^"J^ ^ and 
(j)n(z) — , + (z) and 4>-{z) can be obtained 

by the method of imaginary-time propagation where the 
time variable r in Eq. ([2]) is replaced with an imaginary 
time. 

The initial state with a given initial imbalance at t = 
is constructed as follow 

^(z;0) = e* Ae ^V^MMz) + V^(0)^fl(«), (7) 

where An(0) = n£,(0) — rifl(O) and A#(0) represent the 
initial population imbalance and phase difference in the 
double-well. The dynamical evolution is calculated by 
solving Eq. @ through real-time propagation method. 

Once ip(z;r) is obtained, the atom fraction in the 
left and right well could be computed by 71l(t) — 

/-oo dz l^( 2; ; r )l 2 and "fl( T ) = /o°° dz l^( z ' T )l 2 - Tne Pop- 
ulation imbalance is An(r) = til(t) — tir(t). The phase 
of condensates in the left well is defined as 



J° oo dzIm[V'(2;r)]p(z;T) 
Vl(t) = arctan— . 

J_ oo dzRe[?/'(z;r)]p(z;T) 



and in the right well 



9r(t) = arctan 



r o °° dzlm[?/>(z; r)]p(z; r) 
r o °° dzRe[*0(z; r)}p(z; r) 



(8) 



(9) 



where p(z;r) = i/j*(z;t)iIj(z]t). The phase difference is 
A6(t)=9 l (t)-6 r (t). 

Moreover, we could get the phase at a certain point at a 
certain time which is different from that of the two-mode 
model. 



(z; t) = arctan 



Im 



z: 



;r)] 



Reefer)] 



(10) 



III. THE PHASE SPACE DIAGRAM 

In following calculations, we set j3 = 10, vt = 5 and 
the dimensionless trap geometry are chosen as vm(z) — 
z 2 + 10cos 2 (7rz/5.2), where q' = 5.2 is calculated from the 
experiment parameter in Ref. (l8| . qo ~ 5/xm. However, 
the dimensionless interacting strength giu is a tunable 
parameter. It varies with particle number of the system, 
N. The time step is 5t — 0.001 and space step 5z = 0.01 
in the calculation. 

The evolution of the population imbalance An and the 
phase difference AO can be calculated from Eq. . Fig- 
ure Q] shows the obtained results for gm — 0.01, with 
the initial states being (An = 0.9 A6 = 0), (0.9, w) 
and (0.5, it), respectively. The solid lines in Fig. [TJa) 
and Hlb) depict evolutions from the state (An = 0.9 
A9 = 0). Both the population imbalance and the phase 
difference oscillate around zero (An = and Ad = 0) 
with time. This denotes the Josephson oscillation mode. 
In the other two cases, An is always positive which means 
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FIG. 1: (Color online) The An-r plot (a) and A9-r plot (b). 
The solid lines, short-dashed lines, and short-dotted lines rep- 
resent the evolutions from the initial state (An = 0.9 A# = 0), 
(An = 0.9 AO = 7r) and (An = 0.5 A8 = tv), respectively. 



that the particle number in the left well is always larger 
that in the right well. This is referred to as the so- 
called macroscopic quantum self-trapping phenomenon 
Q. As shown in Fig. [ljb), the short dashed- line oscil- 
lates around A6 = it while the short-dotted line indicates 
that the phase changes monotonously with time. Such 
tendency characterize the difference between the 7r-mode 
MQST and running phase MQST. 

By combining evolutions of An and A9, the An-A9 
phase space diagram can be constructed. Fig. [2ja) dis- 
plays the three kinds of evolutions discussed above, (i) 
The Josephson oscillation mode characterized by the 
time-average particle imbalance (An) r = is represented 
by closed orbits around A9 = (or A9 = 2kir where k is 
an integer), (ii) The 7r-mode MQST is represented by the 
closed orbits around A9 = tt (or A9 = (2k ± 1)tt). (iii) 
The open orbits between the above two types of closed 
orbits refer to the running-phase MQST. The particle 
imbalance (An) T ^ for both the two MQST modes. 
The An-A9 phase space diagram has been produced 
reviously based on the two-mode model As Ref. 
suggested, each orbit denotes a constant energy con- 
tour line. The 7r-mode has relative high energy. The 
energy maxima lie at the center of the 7r-mode closed or- 
bits and there are two maxima symmetrically located on 
both sides of An = axis. On the contrary, the Joseph- 
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FIG. 2: (Color online) The An-A8 phase space diagram for different interaction strengths, (a) gio = 0.01, (b) gm — 0.005, 
(c) giD = 0.001. The orbits are presented by constant energy lines. 



son oscillation mode has smaller energy and the energy 
minimum is located at the points (An = 0, AO = 2fc7r). 

Nevertheless, it is suggested that the two-mode model 
has to be improved to agree with experimental observa- 
tions [18|. Meanwhile, it is shown that the TDGP equa- 
tion provides more accurate simulation of dynamical evo- 
lution of the double-well condensate than the two-mode 
model [iH [23| . Unlike the two- mode model which ig- 
nores the spatial information, the TDGP equation deals 
with the spatial distribution of density and phase, thus 
it achieves quantitative accuracy of dynamics of the con- 
densate. The present study produces the An-AO phase 
space diagram systematically within the framework of 
TDGP equation. 

MQST is a kind of nonlinear effect which arises because 
of the interaction between bosons and might disappear 
if the interaction becomes weaker. In Fig. [2jb) where 
9id = 0.005, the running-phase mode disappears and 
the 7r-mode MQST is surrounded by large closed circles 
referred as the 7r-mode Josephson oscillation. In Fig. (He) 
where g\o = 0.001, the 7r-mode MQST disappears also. 

AO is the conjugate variable of An. We can calculate 
the average phase in each well according to Eqs. ((H]) 
and ©. Fig. [3] shows Ol, 9r and AO for a Josephson 
oscillation case. It is worth noting that the phase of the 
condensate in each well varies much faster than the phase 
difference between the two wells. The phase in each well 
changes almost synchronously with slightly different high 
frequencies. In the two-mode model, 0l and Or are not 
defined. 



IV. THE DISSIPATION EFFECT 

To proceed, we investigate dynamical behaviors of the 
double-well condensate subject to unavoidable dissipa- 
tion. We consider a system described by phase space 
diagram in Fig. HJa) and suppose that it is initially in a 
7r-mode MQST state, (An = 0.9, AO = vr). 

Fig. [4] shows the evolution trajectories calculated us- 
ing Eq. (U). The dissipation parameter is chosen to be 
7 = 0.01, 0.03 and 0.1, respectively. In all the three cases, 
the trajectory passes through three different regimes. In 
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FIG. 3: Average phases in the left well, 8l (a), the right well, 
8r (b) as well as the phase difference A6 (c) for a Joseph- 
son oscillation case corresponding to the long-dashed line in 
Fig. [TJi. 



the 7r-mode MQST regime, the system moves along an 
outward spiral path. Then it enters the running-phase 
MQST regime where the phase difference changes sig- 
nificantly. At last, it slides into the Josephson oscilla- 
tion regime, along an inner spiral path. As discussed 
above, the energy decreases from the 7r-mode MQST to 
the Josephson oscillation. Therefore such trajectories im- 
ply that energy dissipates gradually. Moreover, the larger 
the dissipation parameter is, the faster the system drops 
into the lowest energy state. 

Marino et al. have extended the two-mode model to 
include the dissipation effect [Hj]. The dissipation term 
they employed is based on the address in Ref. |3[ , where 
it was assumed to be an Ohmic current of normal atoms 
which causes energy dissipation. Present study is quali- 
tatively consistent with results of Ref. [32J . Since our ap- 
proach is an extension of the TDGP equation approach, 
we believe that the obtained results provide a more quan- 
titatively accurate description of the dissipation process 
for the double-well condensate and arc expected to be 
compared with experiments, as discussed in Sect. IIIII 

To illustrate the dissipation process clearly, we calcu- 
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FIG. 4: (Color online) Evolution trajectories of a system initially in the 7r-mode MQST state (An = 0.9, AO = n) with different 
dissipation parameters, 7 = 0.01 (solid line), 7 = 0.03 (short-dotted line), and 7 = 0.1 (short-dashed line), respectively. 
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FIG. 5: (Color online) Evolution of An (a) and A8 (b) start- 
ing from a 7r-mode MQST point (An = 0.9, A9 = n) with 
dissipation parameters 7 = 0.03 (short-dotted lines), com- 
pared with case without dissipation (solid lines). 



abruptly from 7r to — 87r, which is the typical character of 
the running-phase MQST regime. When Ad starts to os- 
cillate around —8tt, the system enters into the Josephson 
oscillation regime. During this stage, both An and AO 
oscillate in a sine-like way, but their amplitudes damp 
with time. Finally, An arrives at the value of zero which 
represents an equilibrium distribution of particles in the 
two wells. 




...50. .. 






T=10 




(c) 


_ 40 






" 20" " 




30 


I 


6 -4 -2 




z 


2 4 ( 



2.5 
2.0 

0.5 
0.0 



700 






(d)iooo 




100' 


T=10 


_ ✓ 


200 
400 







FIG. 6: (Color online) Snapshots of phase distributions in the 
double- well at different given times, (a), (b) and (c) depict 
the Josephson oscillation, running-phase MQST and 7r-mode 
MQST cases without dissipation, respectively, (d) shows the 
dissipation prosess from the 7r-mode MQST state with 7 = 
0.03. 



late time-evolutions of the population imbalance and the 
phase difference, taking the 7 = 0.03 case as an exam- 
ple, as shown in Fig. [5] When t is less than 350, the 
system is still in the 7r-mode MQST regime where both 
An and AO oscillate around the center of the 7r-mode 
MQST region, although their amplitudes grow slightly. 
During the interval from r ~ 350 to r ~ 600, AO drops 



An advantage of the TDGP equation approach is that 
it provides spatial information of the phase distribution 
in each well [2iq . Fig. [5] illustrates the phase distribution, 
0(z; r), for different dynamical processes. In all the cases, 
the phase keeps almost a constant in each well, except at 
the edges of the well. Even if the phase varies with the po- 
sition, the deviation is usually very small. This accounts 
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FIG. 7: (Color online) Snapshots of particle density distribu- 
tions in the double-well, (a) depicts the 7r-mode MQST case 
without dissipation, (b) shows the dissipation prosess from 
the vr-mode MQST state with 7 = 0.03. 

for the validity of the use of two-mode approximation. 
An interesting observation is that the phase distribution 
in the well with more particles seems flatter than that in 
the well with less particles. This point is demonstrated in 
Fig-HD 3 ) andH£c) which depict the MQST cases with the 
left well trapping more particles. We also find that the 
phase distribution in the dissipative case is much flatter, 
as if the dissipation tends to suppress phase fluctuation 
in each well, as shown in Fig. [BJd). 

Fig. [B] also shows that the phase 8(z; r) may change 
abruptly near the barrier, but not right at z = 0. Tak- 
ing Fig. [5(c) which describes the 7r-mode MQST as an 
example, there is an insanely steep gradient near z w 0.5 
on each lines. 

We now look at the particle density distribution. 
Fig. [7(a) displays the 7r-mode MQST case without dissi- 
pation, where n(z) exhibits two peaks located in the two 
wells respectively. Obviously, the left peak keeps higher 
than the right one all the time. We note that the min- 
imum between the two peaks does not lie at the center 
of the barrier. It is biased towards the right side of the 



barrier. This accounts for the reason why 8(z; r) begins 
to change abruptly on right side of the barrier. We judge 
that the abrupt change in 6{z] r) takes place at the min- 
imum point of the particle density distribution between 
the two density peaks. Fig. [7(b) shows the dissipation 
process from the 7r-mode MQST with 7 = 0.03. Then the 
system evolves into the running-phase MQST regime (see 
the r = 400 line) and the Josephson oscillation regime 
(see r = 900, 1000 and 2000 lines). At r = 2000, the 
system becomes almost stable. 



V. CONCLUSION 

In conclusion, we have investigated dynamic proper- 
ties of the double-well Bose-Einstein condensate subject 
to dissipation in the framework of time-dependent Gross- 
Pitaevskii equation. The phase space diagram for the 
system without dissipation are systematically produced 
and three typical dynamical regimes of Josephson oscilla- 
tion, 7r-mode MQST and running phase MQST are shown 
clearly. The energy dissipation processes are studied by 
introducing a phenomenological dissipation term to the 
TDGP equation. The evolution trajectories evolve grad- 
ually from the high energy states to the low energy states, 
suggesting that this approach is applicable to describ- 
ing energy dissipation effect. Relevant results have been 
studied previously based on two-mode model. We ex- 
pect that the TDGP approach provides a more quantita- 
tively accurate description of the dissipation process for 
the double-well condensate and the results are expected 
to be compared with experiments. 

Moreover, the TDGP approach produces spatial in- 
formation of the phase distribution and density profile 
which are not directly reflected in two-mode model cal- 
culations. The phase of condensate varies slowly in space 
in each well and may exhibit abrupt changes near the 
barrier, but not certainly at the barrier center. These 
sudden changes occur at the minimum position in parti- 
cle density profile. We also note that the average phase 
in each well varies much faster with time than the phase 
difference between two wells. 
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